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BACKGROUND OF THE INVENTION 

1 . Field of the Invention: 

The present invention relates to a method for detecting circles, and more particularly a 
10 method for detecting circles in an image using horizontal and vertical scanning. 

2. Description of Prior Art: 

Detecting a circle in an image is a basic task needed for computer vision, for example, 
in Surface Mounted Device (SMD) inspection applications. Typically, a Hough transform is 

15 used for computerized circle and curve detection. The principle of the Hough transform is to 
detect parameterized curves by mapping image edge pixels into manifolds in a parameter 
space and finding peaks using a multidimensional histogram procedure. This method, 
however, is expensive in terms of computation and memory needs. Typically, due to memory 
storage limitations, this method also has discretization error in both the image and the 

20 parameter space. 

Recently, several methods have been proposed to improve Hough transform 
techniques in terms of efficiency, accuracy and memory storage. One modified scheme uses 
gradient orientation from edge points. The probablisitic Hough transform attempts to reduce 
redundant information by sampling image data in various ways. The randomized Hough 

25 transform is proposed to detect a curve with n parameters by randomly picking n pixels and 
mapping them into one point in the parameter space, instead of transforming one pixel into a 
hyperplan in n-dimensional parameter space. Other methods, such as the generalization 
Hough transform and the decomposed Hough transform, have been proposed to improve 
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Hough transform techniques. However, the above-mentioned methods are Hough 
transformation based and are therefore, generally slow. 

Therefore, a need exists for a fast, accurate, and robust method of dominant circle and 
ring detection in a given region without transforms (needing less memory). 

5 

SUMMARY OF THE INVENTION 

In one embodiment of the present invention, a method is disclosed for determining a 
circle in a region of interest. According to the method, a first edge point pair is extracted 
along an x-axis of the region of interest, and a second edge point pair is extracted along ay- 

10 axis of the region of interest. The method determines an intersection of a first and second line 
extending perpendicular from a pair of midpoints of the first and second pair of edge points 
respectively, and determines a radius from the intersection to any edge point. The x-axis and 
the>>-axis intersect within the circle. 

The method defines the region of interest wherein the circle is completely contained 

15 within the region of interest, and the circle is the dominant feature within the region of 
interest. 

In one embodiment of the present invention, the method verifies the edge points by 
scanning the image along the x-axis of the region of interest, scanning the image along the y- 
axis of the region of interest, performing a horizontal gradient and a vertical gradient along 
20 the x and ^y-axis of the region of interest respectively, and determining whether a local 
maximum along the gradients match the coordinates for any edge point. 
The match is determined by searching from each edge of the region of interest inward to 
determine the local maximum, and determining whether the local maximum matches the 
coordinates for any edge point. 
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In another embodiment of the present invention, the method determines the match by 
determining a position, relative to the circle of a foreign structure: if the foreign structure are 
inside, searching from each edge of the region of interest, inward, to determine the local 
maximum; if the foreign structure are outside, searching from the center of the region of 
5 interest, outward, to determine the local maximum; and if the foreign structure is both inside 
and outside the circle, randomly scan at a plurality of points each axis of the region of interest 
within Vi of the total axis length, conduct a gradient operation at each point scanned, 
determine a plurality of potential edge points, calculate a coordinate for center of the circle, 
save the coordinate in an array, determine, based on the plurality of edge points a median 
10 value for each coordinate of the center of the circle. The method further determines whether 
the local maximum matches the coordinates for any edge point. 

The region of interest can be selected manually. 

In still another embodiment of the present invention, a method for determining a 
circle in an image is disclosed wherein the method extracts a first pair of edge points along an 
15 x-axis of the image, extracts a second pair of edge points along aj/-axis of the image, and 
determines an intersection of a first and second line extending perpendicular from a pair of 
midpoints of the first and second pair of edge points respectively. 

The method also determines a radius from the intersection to any edge point. 
According to the present embodiment, the x-axis and thej>-axis intersect within the circle. 
20 According to another embodiment of the present invention, the image contains the whole 
circle, and the circle being the dominant feature within the image. 

The method verifies the edge points by scanning the image along an x-axis, scanning 
the image along a^-axis, performing a horizontal gradient and a vertical gradient along the x 
and>>-axis respectively, and determining whether a local maximum along the gradients match 
25 the coordinates for any edge point. The match is determined by searching from each edge of 
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the image inward to determine the local maximum, and determining whether the local 
maximum matches the coordinates for any edge point. 

In still another embodiment of the present invention the method determines a match 
by determining a position, relative to the circle of a foreign structure: if the foreign structure 
5 are inside, searching from each edge of the image, inward, to determine the local maximum; 
if the foreign structure are outside, searching from the center of the image, outward, to 
determine the local maximum; and if the foreign structure is both inside and outside the 
circle, randomly scan at a plurality of points each axis of the image within 14 of the total axis 
length, conduct a gradient operation at each point scanned, determine a plurality of potential 

10 edge points, calculate a coordinate for a center of the circle, save the coordinate in an array, 
determine, based on the plurality of edge points a median value for each coordinate of the 
center of the circle. The method further determines whether the local maximum matches the 
coordinates for any edge point. 

According to another embodiment of the present invention, a computer program 

15 product is disclosed including a computer usable medium having computer readable program 
code embodied therein for causing the computer to detect a circle within an image, the 
computer readable program code. The computer program product includes computer readable 
program code for causing the computer to extract a first pair of edge points along an x-axis of 
the image, computer readable program code for causing the computer to extract a second pair 

20 of edge points along aj-axis of the image, computer readable program code for causing the 
computer to determine an intersection of a first line and a second line extending 
perpendicular from a midpoint of the first and second edge points respectively, and 

computer readable program code for causing the computer to determine a radius from 
the intersection to any edge point. 
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According to one embodiment, the computer usable medium also includes computer 
readable program code embodied therein for causing the computer to verify the edge points, 
the computer readable program code. The computer program product includes computer 
readable program code for causing the computer to scan the image along an x-axis, computer 
5 readable program code for causing the computer to scan the image along aj>-axis, computer 
readable program code for causing the computer to perform a horizontal gradient and a 
vertical gradient along the x and^y-axis respectively, wherein the x and^y-axis intersect within 
the circle, and computer readable program code for causing the computer to determine 
whether a local maximum along the gradients match the coordinates for any edge point. 

10 According to yet another embodiment of the present invention, the computer readable 

program code for causing the computer to determine a match further includes computer 
readable program code for causing the computer to search from each edge of the image, 
inward to determine the local maximum, and computer readable program code for causing the 
computer to determine whether the local maximum matches the coordinates for any edge 

15 point. 

According to another embodiment, the computer readable program code for causing 
the computer to determine a match further includes computer readable program code for 
causing the computer to determine the position, relative to the circle of a foreign structure: if 
the foreign structures are inside, searching from each edge of the image, inward, to determine 

20 the local maximum; if the foreign structures are outside, searching from the center of the 
image, outward, to determine the local maximum; and if the foreign structure is both inside 
and outside the circle, randomly scan at a plurality of points each axis of the image within X A 
of the total axis length, conduct a gradient operation at each point scanned, determine a 
plurality of potential edge points, calculate a coordinate for a center of the circle, save the 

25 coordinate in an array, determine, based on the plurality of edge points a median value for 
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each coordinate of the center of the circle. The present embodiment also includes computer 
readable program code for causing the computer to determine whether the local maximum 
matches the coordinates for any edge point. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Preferred embodiments of the present invention will be described below in more 
detail, with reference to the accompanying drawings: 

Fig. 1 is a diagram of a region of interest having a circle with parameters according to 
one embodiment of the present invention; 

Fig. 2 is a flow chart of a method of circle detection and verification according to 
another embodiment of the present invention; 

Fig. 3 is an illustration of two circles within images according to yet another 
embodiment of the present invention; 

Fig. 4 is an illustration of two additional circles within images having noise according 
to one embodiment of the present invention; 

Fig. 5 is an illustration of a circle within images having noise according to one 
embodiment of the present invention; 

Fig. 6 is an illustration of four circles within images having noise according to still 
another embodiment of the present invention; 

Fig. 7 is an illustration of an image having two circles within the same image; 

Fig. 8 is an illustration of detected circles within an image according to one 
embodiment of the present invention; and 

Fig. 9 is another illustration of detected circles within an image according to an 
embodiment of the present invention. 




DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

According to the geometric properties of a circle, a line that intersects the center point 
of a segment connecting two edge points on a circle, perpendicular to the segment, will also 
intersect the center of the circle. According to one embodiment of the present invention a 
5 method is disclosed for dominant circle detection utilizing the relationship above. The 

following disclosure will be described with respect to Fig. 1. The circle to be detected 106 is 
scanned horizontally and vertically through the center of an image 108 or a region of interest 
(ROI), 0(/tx, vy). The coordinates for the center of the circle 1 10, (xo, yo), can be obtained 
from the center of horizontal segment and vertical segment as shown in Fig. 1. 
10 In Fig. 1, edge points Pl(xi, h y ), P2(x 2 , h y ) are end points of a horizontal segment and 

the edge points Ql(v x , yi), Q2(v x , y 2 ) are end points of a vertical segment. The following 
results can be obtained: 




Where (xo, yo) are the coordinates for the center 1 10 of the detected circle 106 and r is the 
radius of the detected circle 106. 

According to one embodiment of the present invention, the circle 106 to be detected is 
20 assumed to be within the ROI. Further, the ROI center 108, 0(/rx, vy), is assumed to be 

located inside the circle 106. In some cases it can be difficult to select a ROI to exclude other 
structures located inside or near the circle in real applications. Typically, ROI's are 
rectangular, though other shapes are possible. In this case, a rectangular ROI is used. The 
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ROI is selected to reduce the presence of foreign structures to increase the efficiency and 
robustness of the present invention. The ROI can be selected manually. 

According to the present embodiment, edge points PI and P2 can be extracted by 
scanning the image horizontally at the position of row hy 104 and conducting a horizontal 
gradient operation. Similarly, edge points Ql and Q2 can be extracted by vertical scanning at 
the position of column vx 102, and conducting a vertical gradient operation. Then by using 
Equations (1) through (3), the circle parameters can be obtained. 

According to another embodiment of the present invention, a method for handling 
different background situations is presented, e.g., when the background possesses noise 
and/or other structures inside or outside the circle. In one embodiment of the present 
invention, a method is disclosed for determining edge points. The following description deals 
with each successive stage of the method in detail with reference to Fig. 2. 

A method for fast dominant circle detection is proposed using horizontal and vertical 
scanning through the center of an image or a region of interest (ROI), OQix, vy). The 
coordinates for the center of the circle, (xo, yo), can be obtained from the center of horizontal 
segment and vertical segment. A horizontal and vertical scan 204 is conducted and results in 
a pixel array. A gradient operation is performed along the horizontal and vertical plains (pixel 
array), typically the x andj>-axis of the ROI. An edge point is the local maximum in a 
gradient array. The potential edge points are determined 206 for the circle by searching the 
gradient arrays for local maximums. A connectivity verification procedure 208 is performed 
to find the true edge points of the circle. The method obtains circle radiuses j\ and 72 210, 
calculated from edge points. Given a number of scanning times, k 212, the method randomly 
scans the image within a row and column range and determines a center xo 212. The method 
chooses the median value to be the candidate center xo 212. The method verifies whether the 
two radiuses 71 and 72 are equal 214, thus verifying the edge points. If one of the four edge 
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points is not detected accurately, the other three edge points can still be used to calculate the 
correct circle parameters 216. 

In the following paragraphs horizontal scanning is described in detail. Vertical 
scanning is similar to the horizontal scanning, and involves similar techniques. Given aROI 
5 (202) m x n for circle detection, where m and n are the size of the rows and columns in the 
ROI respectively, a horizontal scan 204 is conducted at the position of row hy and results in a 
pixel array A^(i), where / e [1, w] . Then a gradient operation along the horizontal plain is 
performed to obtain a gradient array Gh y (0, where: 

Ghy(i)=\Ahy(i)- Ahy{i-\)\ (4) 

"fi 10 In general, an edge point, or feature point, is the local maximum in a gradient array. If 

m the background is clean, without any noise, and no structure is located outside the circle, the 
*J3 first local maximum is the left edge point 206 of the circle when searching for the gradient 
^ array from left to right. Similarly, the first local maximum is the right edge point 206 when 
1 ~J we search from right to left. However, if noise, e.g., speckle noise, exists in the background, 

15 the first local maximum may be a noise point. Therefore, a connectivity verification 
m procedure 208 can be performed to find the real edge points of a circle. 

Each edge point of a circle has at least two adjacent edge points, where "adjacent" 
means connected. By checking the neighborhood feature points of the first local maximum, 
isolated noise points, or pairs of noise points, can be eliminated from the candidate edge 
20 points of a circle. For low noise background, a 3 x 3 window can be implemented for 

neighborhood connectivity verification. For high noise background, a5x5or7x7 window 
can be implemented. While the window size does not change, the size of columns and rows 
will decrease as the number of columns and rows increases (assuming the size of each 
column and row is equal to every other column and row respectively). Therefore, by 
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increasing the number of columns and rows, the method decreases the likelihood that noise 
will be found in any one column or row during scanning. 

The pseudo-code of connectivity verification 208 is shown below (for horizontal 
scanning): 

5 Scan at row hy to get a row array Ah y (i) 

Conduct gradient operation to get a gradient array Ghy(i) 

For (i=0; i<n/2; z+ +) //search from left to right to find the left edge point 

{ 

If a local maximum is detected Then 
10 If the connectivity of its neighborhood is TRUE, Then 

Exit //obtain the left edge point 

Else 

Continue 

} 

Ol5 For (i=n-l; i>n/2; i—) //search from right to left to find the right edge point 

iM If a local maximum is detected Then 

4! If the connectivity of its neighborhood is TRUE, Then 

s +: Exit //obtain the right edge point 

fi20 Else 

Continue 

} 
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= b It is also possible that some foreign structures are located inside or outside the circle. 

\}25 If the structures are only located inside or outside the circle, then the method can detect the 
circle by conducting outer to inner scanning or inner to outer scanning, respectively. For 
example, if all foreign structures are located outside the circle, by searching for the local 
maximum from inner circle to outer circle, the first local maximum wit h connectivity ^ 11 be 
an edge point. The pseudo-code for two scanning modes according to one embodiment of the 
30 present invention are shown below (for horizontal scanning only): 
If outer to inner mode, Then 

{ 

For (i=0; i<n/2; i++) 

Find the left edge point 
35 For (i=n-l; i>n/2; i-) 

Find the right edge point 

} 

If inner to outer mode, Then 
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{ 

For(i=n/2; i>0; i~) 

Find the left edge point 
For (i=n/2; i>n; i~) 
5 Find the right edge point 

} 

However, if foreign structures are located inside and outside the circle 

simultaneously, the method performs random scanning around the center of ROI. For 

horizontal scanning, given the number of scanning times k 212, the method randomly scans 

10 the image within the row range of [/iy-5, /iy+5] with 5 < m/4. For each scanning array, the 

method extracts the left and right edge points. The center xo is calculated using Equation (1). 

The center xo is written to a new array X. The method also sorts the array and chooses the 

median value to be the candidate center xo 212. The pseudo-code of random horizontal 

scanning is shown below: 

1 5 For (i=0; i<k; i++) 

{ 

Random select a row value t, which t e [hy - 8 ,hy+ 8] 
Horizontal scan at row t to get a row array At(i) 
Conduct gradient operation to get a gradient array Gt(i) 
20 Find the left edge point and right edge point 

Calculate x 0 
X[i]=x 0 

} 

sort the center array XfiJ 
25 x 0 =X[k/2] 

As stated above, according to Equation (3), the method obtains circle radiuses yi and 

72 2 1 0, calculated from PI, P2 and Ql, Q2 separately and respectively. The method verifies 

whether the two radiuses 71 and 72 are equal 214. If foreign structures, such as long lines, 

exist inside and outside the circle simultaneously, it is possible that the two radiuses 71 and 72 

30 are not equal. This is because one or more edge points may not be real or true. However, if 

only one of the four edge points is not detected accurately, the other three edge points can 

still be used to calculate the correct circle parameters 216. The pseudo-code of circle 

parameter verification is shown below: 
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} Wkite(SUCCESS^FALSE)Or(not allpmms^re^ieckedf 



If(the calculated edge point is a real edge point of circle) Then 



If two or more edge points are not detected correctly, a Hough transform based approach can 
15 be used. 

Compared with the conventional Hough transform and other circle detection methods, 
the present invention has the following advantages: 

High computational speed, typically the method does not conduct Hough transform. 
The computational speed is faster than the conventional Hough transform and other circle 
20 detection approaches. Running on a 200MHz computer, the average execution time of the 
circle detection method is about 15|j,s. 

Less memory, the present method does not conduct a transform, all computations are 
focused on scanning several lines in the ROI, therefore, less memory is needed. 



25 parameters without discretization. 

Ring Detection, since that are two scanning modes in the present invention, it can be 
extended to detect a ring. It can also be used to verify whether two circles have the same 
center. 



30 visual C++/MFC with a friendly user interface. Figs. 3 through 7 show the circle detection 
results of artificial images. Fig. 8 and Fig. 9 show the circle detection results of real images. 



High parameter resolution, the present method uses a method which obtains circle 



The present invention can be implemented in any computer language, for example, 
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According to one embodiment of the present invention/a broken line indicates a circle within 
the image located by the method. Circles may be indicated in other ways, for example, by 
lead lines or shading. 

Fig. 3 a shows an ideal circle image. Using the present invention, the circle can be 
5 detected. Fig. 3b shows a circle with speckle noise. By checking the connectivity of edge 
points, the present invention detects the circle. A circle with foreign structures only outside 
(inside) the circle is shown in Fig. 4a (4b). Using inner to outer and outer to inner mode of 
scanning, the present invention can detect these circles simultaneously. The circle can be 
detected through random scanning horizontally and vertically. In Fig. 6, large structures exist 

10 inside and outside the circle. As a result, one of the edge points is not detected correctly. The 
misdetected edge point in Fig. 6a is the left edge point, in Fig. 6b the top edge point, in Fig. 
6c is the right edge point, and in Fig. 6d is the bottom edge point, respectively. Using the 
other three edge points, the method calculates the circle parameters and detects the circle. 
Fig. 7 shows an example of ring detection using two scanning modes. 

15 A method for fast dominant circle detection is proposed using horizontal and vertical 

scanning. Compared to the conventional Hough transform and its derived methods for circle 
detection, the present invention has several advantages, for example, computational speed. 
The proposed circle detection method is also very robust to deal with noise and foreign 
structures. It is suitable for many industrial inspection applications, such as the SMD and 

20 machine parts inspection, as well as computer vision generally. 

Having described embodiments for a method of detecting a circle within an image, it 
is noted that modifications and variations can be made by persons skilled in the art in light of 
the above teachings. It is therefore to be understood that changes may be made in the 
particular embodiments of the invention disclosed which are within the scope and spirit of the 

25 invention as defined by the appended claims. Having thus described the invention with the 
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details and particularity required by the patent laws, what is claimed and desired protected by 
Letters Patent is set for in the appended claims. 
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